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NOMENCLATURE 

= thermal diffusivity of inflation gas, k / p c p , m 2 /s 
= Stefan-Boltzmann constant, 5.729 x 10 ' 8 W/m 2 K 4 
= emissivity of Ballute 

= thermal conductivity of inflation gas, W/m K 
= convective heat flux , W/m 2 
= inflation gas conduction heat flux, W/m 2 
= temperature of free stream, K 
= membrane temperature, K 
= inflation gas initial temperature, K 
= inflation gas bulk temperature, K 
= heat capacity of inflation gas, Joules/kg K 
= free stream velocity, m/s 
= radius of Ballute ring, m 
= inflation gas density, kg/m 3 
= free stream density, kg/m 3 
= time step, s 
= time step index 
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SUMMARY 

A trailing ballute drag device concept for spacecraft aerocapture is considered. A thermal 
model for calculation of the ballute membrane temperature and the inflation gas tempera- 
ture is developed. An algorithm capturing the most salient features of the concept is imple- 
mented. In conjunction with the thermal model, trajectory calculations for two candidate 
missions, Titan Explorer and Neptune Orbiter missions, are used to estimate the stagnation 
point temperature and the inflation gas temperature. Radiation from both sides of the 
membrane at the stagnation point and conduction to the inflating gas are included. The re- 
sults showed that the radiation from the membrane and to a much lesser extent conduction 
to the inflating gas, are likely to be the controlling heat transfer mechanisms and that the 
increase in gas temperature due to aerodynamic heating is of secondary importance. 


INTRODUCTION 

NASA’s Science Mission Directorate in recent years has supported basic research and 
technology investigations for advanced space propulsion through its In-Space Propulsion 
(ISP) office as part of the Solar System Exploration Program. The goal is to develop tech- 
nologies that support reduced trip time with increased payload capabilities, while reducing 
the costs associated with propulsion systems applicable to earth orbit, planetary and inter- 
stellar deep space applications. The intent is to identify, advance, and mature technologies 
that enable NASA exploration missions to the outer planets. 

Research to develop critical technologies for aerocapture provides the impetus for the pre- 
sent work. Aerocapture is a spacecraft flight maneuver performed upon arrival at a planet. 

It utilizes the atmospheric drag to decelerate the spacecraft for capture by the planet's 
gravity into orbit during a single pass. Aerocapture can significantly reduce mission mass 
and cost in comparison to conventional orbital insertion methods [1], 

The aerocapture maneuver can be accomplished with two basic types of systems. The 
spacecraft can be covered with thermal protection material or the vehicle deploy an aero- 
capture device, such as an inflatable heat shield or an inflatable trailing ballute, a combina- 
tion balloon and parachute made of thin, durable material towed behind the vehicle. The 
former is a variation of the conventional aeroshell, using low mass and low thermal conduc- 
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tivity materials with sufficient thermal and mechanical properties to allow orbital insertion. 
The aeroshell can accommodate aerodynamic control to overcome atmospheric uncertain- 
ties and initial entry conditions, while protecting the payload from hypersonic aerothermal 
heating. The present research is concerned with the latter technique, inflatable decelerator 
technology as an alternative to or in addition to the aeroshell, where a large structure can 
remain deflated and stowed away during most of the mission until it is deployed. The toroi- 
dal ballute concept, the subject of this investigation, envisions ballutes, which are large 
enough to allow the spacecraft shock wave to pass through their open center. A spacecraft 
ballute system is schematically depicted in figure 1 . The primary objective here is under- 
standing and insight into the physics, ensuring the functional simplicity of the algorithm and 
recognizing the essential aspects of the phenomenon significant for wide applications. 

Two design reference missions, the Titan Explorer and the Neptune Orbiter are analyzed in 
the present study. Radiation from a strong shock layer to a blunt body can in general be 
significant when CN gas products are present in high concentrations. Titan’s atmosphere is 
more than 98 percent Nitrogen, but contains minor amounts of Carbon, while Neptune’s 
atmosphere contains no Nitrogen. In this study, the essentially cylindrical local geometry 
is open to free gas flow and is not likely to produce significant amounts of CN or shock 
layer radiation at the stagnation region of the membrane. Calculations associated with gas 
composition, state and spectral radiation exchange are beyond the scope of the present 
study. Radiation from the shock layer is, therefore, assumed to be negligible and the sec- 
ondary shocks effects associated with riser cables are assumed to be insignificant. The 
critical thermal consideration here is the temperature attained at the stagnation point of the 
membrane material. 


Aerocapture Ballutes 

Aerocapture ballutes, a form of Inflatable Aerodynamic Decelerators (IAD), are large area 
structures and lack the maneuvering capabilities of aeroshells, since they can produce rela- 
tively little lift. They are deployed at high altitudes in a planetary body's atmosphere to per- 
form deceleration due to their high drag, and are structurally and thermally less robust than 
aeroshells. Although they have been investigated since the late nineteen sixties [2-4], they 
are in the early stages of technology development. One of the challenges is the survival of 
the lightweight ballute material from a mechanical, as well as, from a thermal point of view. 
Toroidal ballutes have been proposed, which are tethered to the spacecraft, inflated and 
later released after sufficient entry velocity has been lost to enable aeocapture. A toroidal 
system can be designed such that all flow processed by the bow shock of the towing 
spacecraft, ignoring influence of the tethers, passes through the hole in the toroid. Signifi- 
cant progress has been made in manufacturing of ballute materials, aerothermal analysis, 
trajectory control, and aeroelastic modeling. Stowing the ballute and its release are doable, 
but the dynamic stability of the spacecraft-ballute system is a more challenging problem. 
These topics, however, are outside of the scope of the present research, since gaining in- 
sight into the thermal problem is the primary goal here. 

Several considerations must be made in choosing the materials for ballute membrane. 
Clearly the internal pressure for inflating and maintaining a certain degree of structural in- 
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tegrity is a factor. Ballutes can be inflated in vacuum prior to atmospheric interface. The in- 
flation pressure, depending on the structural design of the ballute, must prevent large geo- 
metrical deformations during maximum drag. Critical membrane tensile strength is 
approached when substantial atmospheric and aerodynamic forces are encountered. The 
critical membrane tensile strength, however, is a function of temperature. The membrane 
and internal gas temperatures must, therefore, be taken into account in the deployment se- 
quence. 

A number of commercially available materials have been suggested for various preliminary 
considerations [5-7], They appear to be generally non-metallic and mostly polymeric in na- 
ture, such as Kapton, Polyester and Polyethylene. Some contain multiple components, for 
example, fabrics for tensile strength and coatings to render them impermeable. A common 
theme between most of them is that they use a single layer of thin materials referred to as 
thin films. Some of these materials are capable of maintaining sufficient integrity up to 
temperatures of several hundreds of degrees Kelvin. 


Ballute Thermal Problem 

Clearly large ballutes cannot be successfully designed to contain rigid or massive struc- 
tures. They will be analogous to thin membrane balloons laced with cables and only capa- 
ble of withstanding moderate accelerations and temperatures. The spacecraft and ballute 
must develop a significant amount of drag, while minimizing heat flux. Convective 
aerothermal heating rates are kept to low levels so that heat can be radiated away without 
ablation at the surface. In general an extensive model of an aerocapture ballute would 
likely involve non-axisymmetric computation of the external flow field including the space- 
craft flow field and possible shock layer radiation. Entry conditions exist where radiative 
heating may be an important contributor to the total heat load, for example, entry into the 
earth’s atmosphere at Mars return velocities and flight under low-density non-equilibrium 
conditions during the aerocapture maneuver [8], Non-equilibrium radiation can be present 
even when equilibrium radiation is insignificant. Magin, Caillault, Bourdon and Laux [9] 
have presented radiative heat flux predictions of Huygens probe during entry into Titan’s 
atmosphere based on the solution of the radiative transport equations in thermo-chemical 
non-equilibrium. The radiative heating during the entry is primarily due to the CN system. 
Computation of shock layer radiation, however, is rather involved and is outside the scope 
of the present report. The problem is further complicated by possible coupling of multi di- 
mensional convection heat transfer in the inflation gas due to oscillations in addition to ra- 
diation [10-1 1], Oscillations and free convection due to temperature gradients in the 
inflation gas are assumed to not be present or be negligible due to stable fluid stratification. 

Free molecular and transitional flow regimes appear to constitute favorable conditions for 
this concept, depending on the trajectory and attitude stability that the ballute and the 
spacecraft can tolerate. Continuum flow regime poses more severe demands on the me- 
chanical integrity of the system. Feasibility of the concept under these conditions has not 
yet been established. The aforementioned conditions are external to the ballute, while bal- 
lute surface radiation and the inflating gas heat transfer processes may be considered in- 
ternal and surface controlled. Analyses of the internal heat transfer problem have not been 
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reported in the literature, although heat transfer to the inflating gas and radiation from the 
interior surfaces must be assessed in the advanced stages of the program as part of a high 
fidelity research and development effort. The stagnation point of a toroidal ballute is sub- 
jected to spatially maximum aerodynamic heating and is the most likely region to fail due to 
over temperature conditions. Since toroidal ballutes are symmetric and typically have large 
diameters in comparison to the diameter of the circle of revolution, only a cross section of 
the circle of revolution needs to be considered. 

The Computational Fluid Dynamics (CFD) results of Gnoffo [12] for the external flow field 
and the results of computations of Flail and Le [11] for two reference missions produced the 
set of mission parameters employed here. They are presented in a later section. The two 
missions, as mentioned earlier, are the Titan Explorer and the Neptune Orbiter. These mis- 
sion parameters are used as a basis for the present study. Heat transfer to the membrane 
at the stagnation point and the internal heat transfer along a local line of symmetry through 
the circle of revolution are the focus of the present model. 


THERMAL MODEL 

The present model is an approximation of the transient response at the flow stagnation 
point, employing a two-dimensional methodology, where conduction resistance in the mate- 
rial and convection in the gas are neglected. The transient nature of this problem stems 
from two distinct factors. Firstly, the trajectory considerations impose a varying free stream 
velocity and density and thus varying heat input. Secondly, the finite mass of the inflation 
gas is heated or cooled throughout the time period. 

As shown in figure 1 , an area on the stagnation point of the ballute membrane material is 
considered. This area is extended into a volume element through the inflation gas to the 
opposite side of the ballute ring, parallel to the free stream flow. The ballute ring has a ra- 
dius, r and the free stream flow is at a velocity, Lh. Since the membrane is very thin and 
despite its relatively low thermal conductivity and finite heat capacity, we can assume that it 
responds to the heat input rapidly and radiates on both sides instantaneously at a uniform, 
but transient temperature, T w . 

Due to the differing response times of the inflation gas and the membrane, one may con- 
sider a quasi-steady conduction and radiation heat balance situation. One may choose a 
small time step, At during which the temperature of the inflation gas is almost unchanged. 
For the initial time step, a balance between radiation and aerothermal convection heat in- 
put, q a can be assumed according to the Stefan-Boltzmann law. 


(1) q a = 2 8 a T w **4 


For the present model we approximate the gas heating, q c by considering the gas bulk 
temperature Tb, increasing according to the conduction of heat to a semi-infinite fluid at rest 
during the initial time step At, as well as, the consequent time steps. We assume that the 
gas can be treated to be uniformly at the bulk temperature, T b of the previous time step, 
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initially for each new time step. The heat input to the gas can be calculated from the heat 
conduction equation 


(2) q c =k(T w -T b )/( -a At)** 0.5 


Here T w is the equilibrium temperature prevailing in the Ballute material as a result of con- 
vective heating, less conduction to the gas, calculated in the prior time step and balanced 
by the radiative cooling 


(3) T w = [ ( q a - q c ) / 2 ea]**0.25 


q a can be calculated from the standard equation for stagnation point convective heating on 
a blunt body 


(4) q a =C(p,/r ) ** 0.5 11**3 


The dimensions of C do not correspond to a dimensionless constant, but the dimensions of 
q a are as set forth in the nomenclature. Clearly the trajectory imposes a varying free stream 
velocity and density. For or a given ballute, the heat input is a function of the velocity for 
the particular altitude and density. 


Reference [1 1] derived the values of constant C from the two dimensional cold wall formu- 
lation aerothermal simulations of reference [12-13], The velocity and density are supplied 
by the trajectory calculations. 


C= 1 .92 E - 4 for Titan Explorer 

C= 5.80 E - 5 for Neptune Orbiter 


The overall heat transfer to the gas can be calculated from numerical integration of the in- 
flation gas conduction heat flux, q c over the time steps. The properties and the pressure of 
the inflation gas are currently not specified, however a good candidate may be Helium, due 
to its high thermal conductivity and low density. Helium is often used by as an attitude con- 
trol gas and can be stored at very high pressure in vessels. Helium does not participate in 
the radiative process and the radiation from the far side of the ballute wall is essentially 
negligible due to the fourth power dependence of surface radiation. A model algorithm is 
shown below, using the variables listed in the nomenclature, where At represents a uniform 
time step. In a later section verification and testing of the present model will be discussed 
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t = n At ; 


For n = 0 : 

T b (0) = Ti 

q c (0) = 0.0 

q a (0) = C(p,(0)/r ) ** 0.5 U. (0) ** 3 
T w (0) = [ q a (0) 12 e a] **0.25; 


Note that Ti may be different than T and can be much higher if the ballute needs to be 
warmer than the extremely cold free stream temperature to allow deployment. 

For n > 0 : 

q c (n) = k ( T w (n-1 )- T b (n-1 ) ) / ( - a At )** 0.5 

q a (n) =C(p.(n)/r ) ** 0.5 U. (n) ** 3 
T w (n) = [ ( q a (n) - q c (n) ) / 2 sa]**0.25 
T b (n) = T b (n-1) + At q c (n)/ 2 r p c p 


A number of other effects may come into play when considering the real world manifesta- 
tion of this concept. The assumption of temperature uniformity in the membrane can be 
removed by considering the heat conduction and heat capacity of the material. However, 
removal of this idealization only marginally improves the accuracy of the results, as com- 
pared with other factors. The ballute wind ward heating drops off with increasing angle 
from the stagnation point around the ballute ring, but it has no effect on the results due to 
symmetry of the problem on the plane of symmetry. That is, derivatives in the lateral direc- 
tion and therefore lateral conduction vanish as the stagnation point is approached. Alterna- 
tively the same is true due to the relative thinness of the membrane compared to its 
essentially infinite lateral dimensions. Lateral conduction in the inflation gas also does not 
come into play, when the stagnation point is considered. However, more critical issue may 
come into play, if significant interactions between shock waves from the tow cables with the 
boundary layers were present. 

As pointed out earlier, in cases where the heat flow to the inflation gas is negligible for all 
time steps, a model balancing the aerodynamic heating with surface radiation is achieved. 
The surface temperature can then be simply calculated from the following equation. 


T w (n) = [C (p, (n) / r ) ** 0.5 U (n) ** 3 12 so]** 0.25 
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TRAJECTORY FOR AEROCAPTURE 

The trajectories for aerocapture of the missions for the Titan Explorer and Neptune Orbiter 
were calculated utilizing the Traj code [15]. Traj is a computer program used for design of 
TPS for direct atmospheric entry and aerocapture trajectories using a three degrees of 
freedom trajectory simulation program. Stagnation point heat transfer can also be calcu- 
lated in Traj, but only for nose bodies with insulated walls, as opposed to the body section 
with un-insulated walls of the present geometry. 

It is necessary to consider actual mission scenarios to identify realistic atmospheric entry 
conditions. By selecting combinations of departure and entry conditions and estimating the 
mission total weight requirements, representative or baseline missions, can be considered. 
As indicated earlier, reference [11] considers several candidate missions and includes a set 
of mission parameters to Titan and Neptune, the Titan Explorer and the Neptune Orbiter. 
The values for the vehicle mass, initial velocity, altitude and entry angle, the ballistic coeffi- 
cient and diameter of the circle of revolution of ballute ring are tabulated below: 


Vehicle mass, kg 
Initial Flight Angle, degrees 
Initial Vehicle speed, km / s 
Initial Vehicle Altitude, km 
Ballistic Coefficient, kg / m2 
Circle of Revolution Diameter, 


Explorer 

Neptune Orbiter 

393 

335 

-34.3 

-10.8 

9.6 

28.7 

1,200 

1,200 

0.63 

0.36 

13 

25 


Achieving aerocapture for the present mission parameters is not a trivial mater, although 
areocapture conditions do exist when precise mission parameters are calculated and are 
truly reached. Trajectory calculations based on these values were carried out using the Traj 
code [15]. For the Titan Explorer an entry angle of -34.3 degrees and for the Neptune Or- 
biter an entry angle of -10.8 degrees produced aerocapture. The velocity data as a func- 
tion of time from the initial conditions are presented in figure 2. The trajectories are 
sensitive to the values of the entry angle for the initial conditions imposed as noted by ref- 
erence [1 1], Traj exhibited a narrower aerocapture entry corridor than was reported in ref- 
erence [1 1], In both cases, high values for the entry angle produce impact conditions, 
while low values produce skip off conditions. The sensitivities are of considerable impor- 
tance to the orbital dynamics problem, which is outside the scope of the present research. 
They also alter the heat transfer, which make very precise comparisons unwarranted, un- 
less extensive future studies are conducted. 
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HEAT TRANSFER RESULTS 

The free stream velocity, temperature and density as functions of time were extracted at 
uniform time intervals from the trajectory solutions and were used as input data in the 
thermal model algorithm described earlier. The initial free stream atmospheric tempera- 
tures for the Titan Explorer and the Neptune Orbiter missions were 175 K and 148 K re- 
spectively. The values for the emissivity for the membrane were set at 0.9, representative 
for the material under consideration, while the inflation gas heat capacity, thermal conduc- 
tivity and thermal diffusivity were 1 .4E4 J/kg-K, 0.2 W/m-K and 1 .55E-04 m2/s respectively. 
Owing to the relatively large volume of the inflation gas, its pressure and temperature are 
expected to remain only slightly above atmospheric values. 

The model calculated the convective, conductive and radiative heat transfer for the stagna- 
tion point and the inflation gas. The convective heat transfer results are presented in figure 
3. As showed therein, the maximum convective heating occurs at 180 seconds and at 230 
seconds for the Titan Explorer and the Neptune Orbiter missions respectively. The mem- 
brane temperature at the stagnation point is presented in figure 4. The maximum mem- 
brane temperature occurs at the above times and reaches values of 780 K and 839 K for 
the Titan Explorer and the Neptune Orbiter missions respectively. The bulk temperature of 
the inflation gas, as shown in figure 5, increased only of order 20 K for both missions, indi- 
cating that the initial gas temperature is more important in determining the final gas tem- 
perature than conduction during reentry. As indicated earlier, these temperatures may be 
too low for proper deployment of stowed ballutes. If higher temperatures are needed, the 
algorithm allows these values to be set to higher temperatures. 

In a simpler model, reference [1 1] assumed a balance between radiation from one side of 
the membrane and aerothermal convection according to the Stefan-Boltzmann law. In that 
study conduction to the inflation gas and radiation from the inner surface of the membrane 
were ignored. The results presented in reference [11] indicate that the maximum mem- 
brane temperature for the Titan Explorer and the Neptune Orbiter missions occur at ap- 
proximately 170 and 200 seconds and reach mean values of 1052 K and 1327 K 
respectively. These maximum temperatures are reached approximately 10 and 30 seconds 
earlier than those predicted in the present study for the Titan Explorer and the Neptune Or- 
biter missions respectively. 

As noted earlier, there are differences between the trajectory sensitivities and therefore the 
actual trajectories for the present report and those used in reference [11], Although aero- 
capture is achieved in both studies for similar initial entry conditions, since the density and 
velocity conditions encountered have sharp gradients, the aero heating and heat transfer 
are also sensitive to the exact conditions of the trajectory. The membrane temperatures in 
both studies are in reasonable agreement, but they cannot be compared precisely unless a 
more extensive study of the trajectory problem is conducted. The present study, however, 
was extended by computations with analogous physical modeling to reference [11], i.e. 
without conduction through the membrane and having radiation only on one side. The wall 
temperatures achieved, also presented in figure 4, for Titan Explorer and the Neptune Or- 
biter missions occur again at approximately 180 and 230 seconds and reach mean values 
of 937 K and 1007 K respectively, showing a favorable trend in approaching the results ref- 
erence [11], 
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VERIFICATION 


It should be noted that the computational algorithm presented here was tested to verify its 
validity. Exact or approximate analytical solutions for the problem at hand do not exist, 
however there are solutions available in the heat transfer literature for some limiting cases, 
which are both physical and are of interest for application to spacecraft. The first case 
considered is for insulated wall, i.e. where the conductivity of the inflation gas is negligible 
or set to zero. This case approximates a non-ablative, but highly insulating thermal protec- 
tion material used in a reentry environment, where the aerothermal heating is balanced by 
surface thermal radiation. The surface temperatures obtained from the present algorithm 
for this case very closely matched the values obtained using the Stefan-Boltzmann law. 

The second test case considered is for constant values of the free stream velocity and 
density, which imposes constant convective heating at the surface. Conduction to the infla- 
tion gas takes place, but the heat transfer processes here are controlled by the highly non- 
linear surface radiation transfer mechanism. In this case, consistent with the missions 
considered for the present analyses, the conductive portion of the total heat transfer from 
the surface was only a few percent, the rest being dissipated by radiation. The free stream 
velocities and densities were set at arbitrary, but physically realistic middle range values for 
the tests cases. Although heat transfer literature contains a large number of rigorous ana- 
lytical solutions involving heat conduction in solids, it is devoid of a suitable or similar case 
for comparison. The very low heat conduction to the inflation gas for the ballute system, 
however, allows an approximate comparison of wall temperatures for a purely radiative sur- 
face to the values obtained from the present algorithm. For this test case, the surface tem- 
peratures obtained from the present algorithm were lower than the values obtained using 
the Stefan-Boltzmann law by less than three percent, a reasonably expected outcome. 


CONCLUSIONS 

The present model predicts the approximate peak ballute stagnation point temperatures at 
roughly the same times as predicted in reference [11], but at somewhat lower values. 
Aerocapture is achieved in both studies for similar initial entry conditions, but since the 
density and velocity conditions encountered have sharp gradients, the aero heating and 
heat transfer are also sensitive to the exact conditions of the trajectory. The wall tempera- 
tures in both studies are in reasonable agreement, however, a precise comparison cannot 
be made unless more extensive studies of the trajectory problem is conducted. The lower 
temperatures may be a result of the differences in the trajectories and the more realistic 
two sided radiation boundary conditions set forth in the present study. The temperatures 
obtain are within the tolerable limits for the membrane materials considered. So far as the 
internal and surface controlled thermal modeling is concerned, the likely oscillatory nature 
of the dynamical spacecraft-ballute system poses a challenging mechanical design. The 
heat transfer is primarily controlled by radiation, although forced convective currents within 
the inflation gas may be generated, possibly lowering the temperature of the heated area. 

In an extension to the present model and depending on the nature and magnitude of the 
flow field outside and inside the ballute, a quasi-steady aerothermal convective and radia- 
tive problem may be constructed, for example for analyses of the Hypersonic or Supersonic 
Inflatable Decelerators HIAD and SIAD, respectively. In such a model the heat conduction 
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resistance in the membrane and heat capacity of the membrane material, although of sec- 
ondary importance, may also be taken into account. These extensions, however, would 
render the present method more complex, requiring further development. 
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Schematic representation of the spacecraft-ballute system 
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Free-stream velocity as a function of time from the initial conditions 



Figure 3. Convective heat transfer as a function of time 
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Figure 4. Membrane temperature at the stagnation point 



Figure 5. Bulk temperature of the inflation gas 
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